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ABSTRACT 

We investigate the possibilities that pulsars act as the lens in gravitational microlensing 
events towards the galactic bulge or a spiral arm. Our estimation is based on expectant 
survey and observations of FAST (Five hundred meter Aperture Spherical Telescope) 
and SKA (Square Kilometer Array), and two different models of pulsar distribution are 
used. We find that the lensing rate is ^ 1 event/decade, being high enough to search 
the real events. Therefore, the microlensing observations focusing on pulsars identified 
by FAST or SKA in the future are meaningful. As an independent determination 
of pulsar mass, a future detection of microlensing pulsars should be significant in 
the history of studying pulsars, especially in constraining the state of matter (either 
hadronic or quark matter) at supra-nuclear densities. The observations of such events 
by using advanced optical facilities (e.g., the James Webb Space Telescope and the 
Thirty Meter Telescope) in future are highly suggested. 
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1 INTRODUCTION 

Pulsar s could be norm al neutron stars or quark 

stars (jLattimer fc PrakashI l2004h . and it should be very 
important to affirm or negate the existence of either neutron 
or quark stars in order to guide physicists in studying the 
nature of fundamental strong interaction. With regard to 
the possible ways of identifying quark stars (e.g.. IXull2008l . 
for a review), it could be very straightforward and clear if 
we would find low-mass quark stars, since neutron stars are 
essentially gravitation-bound while low-mass quark stars 
are mainly confined by strong interaction. If we can detect 
a pulsar-like star with mass < O.IM©, then it is surely a 
quark star. 

How to measure the mass of stars? Up to now only 
the masses of compact stars in binary systems have been 
determined, by either Keplerian or post-Keplerian parame- 
ters. The lowest mass detected of an eclipsing X-ray pul- 
sar, SMC X-1, is 1.06tg;}J Mq, which is near the mini- 
mum mass expected for a neutr on star produced in a super- 
nova (|van der Meer et al.ll2007l ). However, most of pulsars 
are isolated, and it is still a big challenge to measure their 
masses. It was suggested to observationally determine an iso- 
lated neutron star's mass from the red-shift (as a function of 
the ratio of mass to radius, M/R) and pressure broadening 
(as a function of M/R^) of an absorption spectrum, yet no 
atomic line has been detected with certainty in the thermal 
X-ray spectra (e.g.. [Xu.200 2'). 

Interestingly, gravitational microlensing (e.g., iMad 



l2008h . which makes use of the temporal brightening of a 
background star due to intervening object, would provide 
us a powerful method to measure the masses of compact 
objects. The amplification A of a background star by the 
passage of a pulsar has a relation of 



A 



where u = t'/Re is a parameter, with th e undeflected dis - 
tance, r', and the Einstein radius, Re (jPaczynskil llOsd ) ■ 
which is almost a function of pulsar mass M and r if 
the background star's distance D ^ r [see Eq.([T|]. One 
can then directly obtain the mass-distribution of the lens 
objects via the light curve (|Glicensteinl booal ). With the 
method of microlensing, we can not only determine the 
mass distribution of ma ssive compact halo objects (MA- 
CHO, |Alcoc^^eF^0^92)_aBdn^^ objects in the Galac- 
tic center region i Wex et al. Ill99^ ). but may also measure 
the masses of isolated pulsars. 

The possibilities of microlensing neutron star have 
been discussed previously and several estimation of lens- 
ing rates had been presente d in the literatures (|Horvathl 
1 19961 : ISchwarz fc Seidell |2002| ). In this paper we reestimate 
the lensing rate by calculating the solid angle Einstein ring 
swept due to pulsar's proper motion, using two different 
pulsar distribution models. In addition, o ur estimation also 
takes exp ectant observations o f FAST (|Nan et al. 1 12006| ) 
and SKA (|johnston et al. Il2007h into account, since the two 
new telescopes, FAST and SKA, will greatly enhance our 
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ability of searching and observing pulsars, and the number 
of pulsars with measured proper motion in the Galaxy is 
crucial to estimate the lensing rate. Lorimer et al. (2006) 
used the results from recent surveys with the Parkes Multi- 
beam system to derive a potentially detectable population 
of 30000 normal pulsars. By using this number and assum- 
ing 3000 potentially dete ctable millisecond pulsars in the 
Galaxy l|Lvne et al. 1 1 199 ^!), Smits et al. (2009a) estimated 
that the all-sky survey with only the 1-km core of the SKA 
located in the southern hemispheric would detect about 
14000 normal pulsars and about 6000 millisecond pulsars. 
Under the same assumption, Smits et al. (2009b) noted that 
more than 7000 previously unknown pulsars could be dis- 
covered by FAST in the Galactic plane (|fe| < 10°) with an 
observation time of 1800 s per-beam. Further more, by regu- 
lar timing these pulsars using FAST and SKA, the accurate 
positions and proper motion of them can also be determined. 
We then propose to search possible microlensing events due 
to pulsars identified by FAST and SKA and to carry out 
microlensing observations coupled with radio observations 
by means of future optical projects (e.g., the James Webb 
Space Telescop«0 or the Thirty Meter Telescope Projec10). 

The lensing rate of our result indicates a high possibility 
of observing microlensing event due to pulsars in the future 
with the help of FAST and SKA. We expect an independent 
measurement of pulsar mass through these microlensing ob- 
servations, especially to discovery low-mass pulsars, in order 
to understand the real nature of pulsars. 
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where D ^ r. We assume that the number density of stars 
in the galactic bulge or a spiral arm is constant in following 
estimation. Because the velocity distribution of pulsars is 
not a function of r and D, one can integrate the velocity 
distribution independently of position distribution. We can 
thus just apply a mean velocity in estimating the lensing 
rate. Hansen and Phinney (1997) suggested a Maxwellian 
distribution for the kick velocities of pulsars. 
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where a — 190 km/s and mean kick velocity Ukick ~ 
300 km/s are chosen. Considering that only the part of 
the velocity vector that lies in the lens plane is effective, 
one s hould use the mean ve locity projection on the lens 
plane (jSchwarz fc Seidelll2002l ). 
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where the projected velocity, v, is order of 200 km/s. 

We use two different models of pulsar number density 
from Hartman et al. (1997), which resemble the models from 
Narayan (1987) and Johnston (1994), see also Schwarz and 
Seidel (2002), in our following simulations. The first model 
of pulsar distribution is 



2 THE LENSING RATE 

We consider a pulsar at distance r with velocity v, microlens- 
ing a background star at distance D. The Einstein ring of 
this pulsar sweeps a solid angle Sn on the celestial sphere 
during a period time of t. The Einstein radius is 
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where _Rs = 2GM/(? denotes the Schwarzschild radius. This 
solid angle Sn depends on the mass M, the distance r to the 
pulsar and also the distance D to the star, 



S^{M,v,r,D) = ^-. 

r r 



(2) 



We expect that a microlensing event occurs when a star 
in the background falls into this solid angle. We assume that 
our telescopes can identify a total of pulsars, and S square 
degrees of the Milky Way is visible to it, then the lensing 
rate per unit time can be estimated as 



P ■ 



N 



(3) 



where A'^star denotes the total number of stars visible. 

If we take pulsar distribution into account, the sum 
turns into multiple integral on the space of position (r and 
D) and velocity (v). Then the lensing rate p can be rewritten 



(7) 



where R is the radial distance of the pulsar to the galactic 
center in the galactic plane, and Rw = 5 kpc. The second 
model is 



, Cp2 , (-R— -Rmax)^ 



(8) 



where Rw = 1.8 kpc and iimax ~ 3.5 kpc. The normaliza- 
tion constant is cp2 — 0.204 for the given choice of -Rmax. 
For the z-dependence we apply 
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with a — 0.45 kpc ()Lvne et al. 1 1l998l : ISchwarz fc Seidell 
|2002| ). 

We transform the pulsar distribution into the spherical 
coordinates (r, 6, (j>) with the sun at the origin via 



r sin 9, 



R^ = cos^ e + RIc - 2ri?sc cos 6 cos < 



(10) 
(11) 



where R is the radial distance of the pulsar to the galactic 
center in the galactic plane and z is the height, and i?sc = 
8.5 kpc. For pulsar distribution npi, we have 
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^ http://www.jwst.nasa.gov/ 
^ http://www.tmt.org/ 
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Figure 1. The lensing rate density as function of 6 and <j> in 
the model of pulsar distribution described by Eq. jTjl, where 6 is 
the azimuth angle and ip is the polar angle. Color bar stands for 
lensing rate in units of number of events/decade/degree ^. 
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where y = D/r. If we define x — r/Rsc, then the expression 
can be reduced to 
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We consider that FAST is expected to detect about 7770 
pulsars in about 70 x 10 square deg rees of the Milky Way of 
observing time in less than a year (|Nan et al. II2006I ). while 
SKA is expected to detect about 15000 pulsars in 290 x 10 
square degrees of the Milky Way. If one chooses AI — IMq, 
r G [0,5kpc], t = 1 yr, iv'= 15000, TV^tar = lO", v = 200 
km/s, and S = 2000 square degrees, then we have p ~ 12 
event/year from Eq. (|13|l . 

Certainly the total number of stars visible is less than 
10^^ if the effect of extinction is included. Based on the pho- 
tometric maps of the galactic bulge released by the OGLE 
(Optical Gravitational Lensing Experiment) project, which 
contain photometry of about 30 million stars from 49 fields 
coveri ng 11 square degrees i n different regions of the galactic 
bulge l|Udalski et al. II2002I ). we estimate that in the region 
of S = 2000 square degrees, 10^ stars would be visible, then 
the lensing rate should be p ^ 1 event/decade according to 
Eq. (|13|l . The lensing rate density as function of 6 and is 
shown in Fig.l. 

For pulsar distribution np2 in the model of pulsar dis- 
tribution described by Eq.®, we have accordingly 
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Figure 2. Same as in Fig. 1, but in the model of pulsar distribu- 
tion described by Eq.||8ll. 
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From the same process as for npi of Eq.(l7]), then we have 
p ^ 2 event/decade if one chooses M = IMq, r G [0, 5kpc], 
t = 1 yr, iV = 15000, A^star = 10", v = 200 km/s, S = 
2000 square degrees, and takes the effect of extinction into 
account. The corresponding lensing rate density as function 
of 9 and is shown in Fig. 2. 

The lensing rates estimated above are for the microlens- 
ing events in the galactic center. Considering that the ob- 
servation of FAST is limited to the spatial arms and that 
pulsars mainly distribute on the disk of the Milky Way, we 
then transform the pulsar distribution into cylindrical coor- 
dinates with sun at the origin via 



R^ = r'^ + Rsc — 2rRsccos(j>. 

Using the cylindrical coordinates, we have 



(15) 
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For pulsar distributions of both npi and np2, we calcu- 
late following the way we did previously and choose Al = 
IMq, r G [0, 5kpc], t = lyr, N = 10000, v = 200 km/s, and 
S = 2000 square degrees. Considering the effect of extinc- 
tion, we choose Astar = 10^. Finally we can also obtain a 
lensing probability of p ^ 1 event/decade for FAST. 



3 A PROPOSAL OF DISCOVERING 
MICROLENSING PULSARS 

We have estimated the probability of observing microlensing 
pulsars by means of FAST and SKA. Besides the encourag- 
ing lensing rate according to our results, a feasible searching 
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strategy is also essential for discovering a microlensing pul- 
sar. We therefore propose to carry out microlensing pulsar 
observation coupled with FAST and SKA observations in the 
future, which are hopeful to determine the mass of isolated 
pulsars. 

Microlensing pulsar observation should be based on the 
expectant data of FAST and SKA. We propose to systemati- 
cally list the information of new pulsars discovered by FAST 
and SKA in the future, including the position, proper motion 
and distance, and high proper motion pulsars should be es- 
pecially focused on. The position of these pulsars in the next 
few decades could also be predicted. Then we could compare 
the predicted position of these pulsars with the position of 
stars in the galactic bulge and spiral arms, and background 
star candidates of microlensing events whose position would 
be quite close to the predicted position of pulsars could be 
picked up. With these candidates, we should formulate a 
plan to monitor certain background stars for the predicted 
microlensing events in a certain period of time in the future. 

Microlensing pulsar observation could be carried out 
on ongoing lensing projects, and future advanced facilities 
are also expected to benefit the observation. The magni- 
fication of a background star could be about 0.3, on the 
condition that a puls ar just arriving at the edge of the Ein- 
stein ring l|Maoll2008l ). and the time scale of the microlensing 
event should be about 15 days. The observation can then 
be done on current telescopes in the infrar ed band where 
the interstellar extinction is much weaker (|Horvath|[l99^ : 
lUdalski et ani2002h . Besides ongoing gravitational lensing 
projects, we could also look forward to future projects, es- 
pecially the Thirty Meter Telescope (TMT) Project. The 
TMT project, which is scheduled for the next decade, will 
greatly promote our study of gravitational lensing on cos- 
mology, galaxy formation and the distribution of lensing 
mass. TMT's great capabilities of resolving smaller, fainter 
source populations will allow a much higher sk y density of 
background sources to be used (|Carlberg||200i ). which will 
directly enhance the possibility of discovering microlensing 
pulsar events. Therefore, TMT together with FAST and 
SKA will provide us a good opportunity to observe mi- 
crolensing pulsar events and determine the mass of isolated 
pulsars. 

We note that the proposed microlensing pulsar observa- 
tion is costless, that is to say, it does not need to monitor a 
target all the time. Our estimation has shown the possibility 
of microlensing a pulsar, and the key of the observation is 
predicting the microlensing events based on the large and 
sensitive database of FAST and SKA, and then monitoring 
microlensing candidates of background stars in a predicted 
short period of time. Even though the prediction can not be 
precise and the lensing rate is relatively small, the costless 
observation is still meaningful and should be carried out. 



4 DISCUSSION AND CONCLUSIONS 

We have estimated the lensing rates of pulsars towards 
the galactic bulge and spiral arms. Our estimation shows 
that, for FAST and SKA, the lensing rates (oc A*') are ;^ 1 
event/decade at least, which are much higher than the esti- 
mations made in previous literatures. The number of pulsars 
with measured proper motion, A'^, is crucial to estimate the 



rate of lensing events. In table 1, we summarize the expec- 
tant number of pulsars to be detected by FAST and SKA, 
based on previous survey simulation results. In addition, the 
lensing rate should increase significant ly if the p opulation of 
rotating radio transients (RRATs, .Keanell2009h is included, 
since the total number of RR ATs is at least sev eral times 
of galactic active radio pulsars (|McLaughlinll2006l 'l. RRATs, 
which can be precisely located by timing and dispersion mea- 
surement, should then be one of the key targets of FAST and 
SKA. 

We also note that the lensing rates we showed are based 
on photometric microlensing events. If we consider astromet- 
ric microlensing, which makes use of the shift of the centroid 
of the combined images of the light source, the cross sec- 
tion would increase by a factor of ~ 10 (|Horvathl [l99^ : 
ISchwarz fc Seidelll2002l ). As our lensing rate is proportional 
to cross section, the astrometric microlensing probability 
would be ~ 10^ times higher than that presented previously 
in §2. It could then be realistic to search astrometric mi- 
crolensing pulsars and to detect the events by facilities with 
very high position precision in the future. 

The lensing rates indicate that it is hopeful to mea- 
sure the mass of an isolated pulsar in the future with the 
method of microlensing. We propose to do catalogue com- 
parison and microlensing prediction for pulsars identified by 
SKA and FAST in the future, and to carry out microlensing 
observation coupled with radio observation in order to de- 
tect microlensing pulsar events and to measure the masses 
of isolated pulsars with advanced optical facilities. 
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Table 1. Summary of previous survey simulation results of FAST and SKA. Numbers in the parentheses represent known pulsars, 
and Dec are the longitude, latitude and declination in the galactic coordinate system, respectively. 





Detectable pulsars 


FAST'^ 


SKA'' 






All Sky 


20° < i < 90° 20° < « < 90° 
|fe| ^ 10° \b\ ^ 10° 


0° < « < 85° & 155° < I < 360° 
\b\ ^ 5° 


Dec < 50° 


Normal pulsar 


~ 30000" 


~ 5700(352) ~ 7000(418) 


~ 11000 


~ 14000 


Millisecond pulsar 


~ 30000'' 


- 550(14) ~ 770(20) 


~ 4000 


- 6000 



Note. - The results are from "a" : iLorime'^ ||2006|'| . "fe": lLvne et al. I l|l998l') . "c": ISmitsl ll2009bl'l . "d": ISmitsl l l2009al') . 
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